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In low-temperature anti-ferromagnetic LaMnOa, strong and localized electronic interactions 
among Mn 3d electrons prevent a satisfactory description from standard local density and gen¬ 
eralized gradient approximations in density functional theory calculations. Here we show that the 
strong on-site electronic interactions are described well only by using direct and exchange correc¬ 
tions to the intra-orbital Coulomb potential. Only DFT+U calculations with explicit exchange 
corrections produce a balanced picture of electronic, magnetic and structural observables in agree¬ 
ment with experiment. To understand the reason, a rewriting of the functional form of the + U 
corrections is presented that leads to a more physical and transparent understanding of the effect of 
these correction terms. The approach highlights the importance of Hund’s coupling (intra-orbital 
exchange) in providing anisotropy across the occupation and energy eigenvalues of the Mn d states. 

This intra-orbital exchange is the key to fully activating the Jahn-Teller distortion, reproducing 
the experimental band gap and stabilizing the correct magnetic ground state in LaMnOa. The best 
parameter values for LaMnOa within the DFT (PBEsol) + U framework are determined to be U = 8 
eV and J = 1.9 eV. 

PACS numbers: 71.15.Mb, 71.20.-b, 75.30.Et, 75.25.Dk 


I. INTRODUCTION 


LaMnCb (LMO) is characteristic of the ABO3 fam¬ 
ily of strongly correlated transition metal oxide per- 
ovskites, which generally exhibit complex phase dia¬ 
grams, as a result of subtle coupling across several 
distinct mechanisms. 1 Bulk, thin film, and interfa¬ 
cial LaMnOs are subject to a multitude of symme¬ 
try breaking mechanisms, including crystal field,® oc¬ 
tahedral distortion,® orbital ordering and Jahn-Teller 
distortion} 3 "^ Mott-type strong d electron Coulomb in¬ 
teractions (direct and exchange) P® and charge transfer 
ordered (Verwey) states® 11 . All of these mechanisms 
are believed to exist and compete in varying ways in this 
material. As a result, LaMnC^ naturally exhibits a rich 
phase diagram as a function of temperature and pressure® 
as well as dopant concentration} 12 ®^which together make 
LaMnC>3 the single most examined metal oxide in the 
LaXC>3 class (where X is a transition metal atom) ! 14 Dop¬ 
ing on the ABO3 A site provides a particularly rich field 
of experimentally observed phenomena, with both Na 
and Ca doped Lai_ a ,A a Mn 0 3 exhibiting colossal mag¬ 
neto resistance (CMR(P^ 16 and a Seebeck coefficient 
that can exhibit positive or negative values which may 
lead to potential thermopower applications!- 17 Pure bulk 
LaMnOs is spin polarized and non-polar, but recent the¬ 
oretical work shows that the magnetic state in Sr doped 
Lai-^Sr^MnOs may be controlled through variation in 
the electric polarization stateP® Recent multi-ferroic the¬ 
ory predicts novel magnetic properties due to t2 g ferro¬ 
magnetic superexchange in Ti doped LMO interfaces.^ 
Finally, the interface between Lai-^Sr^MnOs and a fer¬ 


roelectric shows a polar state that also has a reversible 
orbital polarization!^ 

This interest in LMO from condensed matter and 
materials scientists underscores the value of a reliable 
first principles description based on, for example, den¬ 
sity functional theory (DFT). In particular the magnetic, 
electronic and crystal structure should be accessible si¬ 
multaneously within a low-cost computational frame¬ 
work. Unfortunately previous Hartree-Fock, DFT and 
hybrid-functional examinations of bulk LMO show that 
that obtaining a satisfactory description is not trivial. 1 20 

In this work we show the limitations and successes 
of two different DFT +/7 methods. The Dudarev et al. 
Coulomb correction}? 1 here called I/ e ff, averages out ex¬ 
change effects of the Mn d shell, and we show that it can¬ 
not simultaneously reproduce the bulk band gap, struc¬ 
ture and magnetism. The dedicated anisotropic exchange 
term in the Liechtenstein et al. Coulomb correctional 
here called U\J, dramatically improves the description 
of LMO. The U\J method answers the specific call for a 
practical DFT-based methodology capable of reproduc¬ 
ing t he g ap, structure and magnetism simultaneously in 
LMOP^This is useful as understanding the coupling be¬ 
tween electronic, magnetic and lattice de grees of freedom 
in LMO is a matter of persistent interest! 2 ! 23 ! 24 ! 

Using the U\J method we show the importance of 
Hund’s coupling in LMO. Intra-orbital exchange can en¬ 
ergetically order the orbitals of the Mn 1 2 S 3 e 5 1 ion, which 
in turn strongly affects inter-orbital magnetism and the 
size of the LMO band gap. The Mn e^ 1 occupancy 
polarizatio d 19 ! 25 ! 
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Ie g cr _ fx 2 —y 2 a f3z 2 -r 2 a /i \ 

~ f 2 2 + fo 2 2 ’ 1 ’ 

for the x 2 — y 2 and 3 z 2 — r 2 occupancy eigenvalues (/) 
where a labels spin, is highly sensitive to intra-orbital 
exchange term J in the U\J scheme. By modifying the 
sign and value of 7r es<T , we correct the DFT description of 
Jahn-Teller (JT) distortion, and the electronic and mag¬ 
netic structures of LMO. In addition the U\J calculations 
provide insight into the origin of magnetic, electronic and 
structural ordering in LMO. 


II. METHODOLOGY 

Periodic plane wave density functional theory ( DFT) 
calculations are performed using the VASP softwarepSI 2 ^ 
the local density approximation (LDA PZ81)Perdew and 
Zunger^, and the generalized gradient approximation 
(GGA) in the form of the Perdew-Burke-Erzenhof solids- 
adapted exchange correlation functional (PBEsol) ! 29 ! 30 ^ 
Valence electrons are descri bed u sing the projector aug¬ 
mented wave (PAW) method 31 32 1 with core states (up to 
4d in La, 2 p in Mn, and 2s in O) frozen at their atomic 
reference states. Plane-waves were cutoff above a kinetic 
energy of 520 eV, and a 5 x 4 x 5 k-point mesh of was em¬ 
ployed for the LaMnOa unit cells. All relaxed structures 
fulfill a convergence criterion of less than 0.01 eV/A, 
for both ionic forces and volume-normalized stresses (as 
standard in VASP). 

DFT has known shortcomings in the prediction of the 
electr onic st ructure of materials with localized electronic 
states!^ 33 * 34 ! A typical example are the bands derived 
from Mn d orbitals in LaMnOa: the errors can be cor¬ 
rected to various extents by employing Hubbard-U type 
corrections to account for intra-atom ic C oulomb interac¬ 
tions in the DFT+U approach! 21 l 22 l 34 l The most pop¬ 
ular and simplest Coulomb correction is the “Spheri¬ 
cally Averaged” scheme of Dudarev et al. pH here called 
DFT+t/ e ff, which has only a single effective U parameter, 
U e g. A more sophisticated approach is the “Rotationally 
Invariant” scheme of Lichtenstein and ZaanenpS which 
we label here as DFT+t/| J. Note the simpler Dudarev 
U e s approach was developed after the Liechtenstein U\ J 
approach, and both are fully rotationally invariant. 

Both DFT+U methodologies add Hartree-Fock type 
corrections to the DFT total energy that act on a local 
sub-space of atomic-like orbitals. The DFT+/7 e g total 
energy is given by 


E'DFT+Lfeff = -®DFT H-££(/--./£) (2) 

at i,cr 

where Edft refers to some chosen flavor of electron 
density-based exchange-correlation approximation (LDA 
or GGA in our work). The index at specifies the Mn 


sites where the correction is performed. The eigen- 
occupations fi cr of the electronic on-site density matrix 
are labeled by spin a and index i which represents a lin¬ 
ear combination of angular momentum quantum num¬ 
bers (which in our case ranges over the five magnetic 
quantum numbers m = —2, —1,0,1, 2 for the 3 d Mn or¬ 
bitals) . U e s = U — J is the Hubbard-type energy param¬ 
eter for this approach while U and J are the separate 
direct and exchange Coulomb parameters* 2 ^ (see also Ap¬ 
pendix A). 

For our work here, the DFT+[/| J total energy is best 
rewritten as an added correction to the DFT+Deff ap¬ 
proach (as detailed in Appendix A) given by 


Edft+u\.j — Adft + E corr 

= -EdfT + £ £(-/W - fio) + 

at i,cr 


r, £] @ij fitrfja' 

<j<j' ,ij 


A X^f ia f ja S a<Tl . (3) 


The correction to the DFT band energy eigenvalue 
steins from the occupancy derivative of the correction 
terms given by 


U<sK ^ ~ 

+ £ C ft 7 fj a < — A X^fjaSaa' , 

3,o' 

where the first term is the U e g correction and the second 
and third terms are the added contribution from the U\J 
scheme. For compactness and for use below, it is useful to 
collect all occupancies or energy eigenvalues for the same 
spin into a vector f a or e a in order write these corrections 
in matrix notation. For atomic d shells, the Appendix A 
shows that 


a corr _ d-^corr 

~ df ta 


Ae“ OTr = t/ eff Q - + JA°f a + JB°f s (4) 

where a represents the opposite spin to a. For canonical 
t 2 g and e g orbitals, the dimensionless matrices A a and 
B a are 


A a = 


and 


B a = 
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Both DFT+ U methodologies permit the description of 
electron localisation phenomena, that stem from Hartree- 
Fock physics and the related removal of self-interaction 
errors, which enable essential long-range ordering (or¬ 
bital, spin, charge and lattice degrees of freedom) P2 For 
Mn in LMO, delocalised s and p orbitals typify the 
weakly correlated electronic states successfully described 
by DFT, while the localised Mn d states require the + U 
correction. In the DFT+£/ e ff approach, U e g in Equation 
([2]) provides occupation-dependent corrections to DFT, 
while the DFT+U] J approach in Equation © adds fur¬ 
ther degrees of explicit spatial/orbital dependent correc¬ 
tions. Both corrections provide a basis for for energy 
splitting of d orbitals (and related symmetry breaking 
and orbital polarization) on top of splittings due to spin 
exchange and/or crystalline geometrical distortions al¬ 
ready present at the LDA or GGA density functional 
level. 

The U\J correction variety in Equation Q is most rele¬ 
vant to materials with strongly interacting electrons with 
an explicit orbital symmetry dependence,^ for exam¬ 
ple, Fe-based superconductors,^ heavy fermion metals 
non-collinear magnetic materials) 38 * 39 * and orbitally or¬ 
dered materials in which Hund’s coupling is critical to es¬ 
tablishing the correct insulating or metallic character.® 
Although the anisotropic exchange corrections to DFT 
have successfully been used to describe manganese oxides 
in the past) 39 * 41 !® we believe our work is the first explicit 
calculation and analysis of the U\J exchange matrix ele¬ 
ments and anisotropic splitting for LaMnO. 3 . 

In our work, our global coordinate system is chosen to 
align the orthogonal x',y,' z' axes along the LMO unit 
cell ( a,b,c ) vectors. A local x,y,z basis for each Mn is 
defined by aligning the local axes with the Mn-0 bonds 
of each tilted MnOg octahedron (see Figures [l] and [4|: 
the local x axis is chosen along the shortest Mn-0 bond 
(strongly JT active), the local y axis along the interme¬ 
diate length Mn-0 bond (here called apical), and the 
local 2 axis is along the longest Mn-0 bond (strongly JT 
active). Use of this local basis is more convenient for 
analysis of the electronic states and occupancies. The 
transformation from global to local coordinates is per¬ 
formed for each relaxed geometry by employing a direct 
polynomial-based transformation of orbitals (detailed in 
Appendix B). Unless specifically noted, orbitals and oc¬ 
cupancies refer to the local basis. 


III. RESULTS AND DISCUSSION 

At 750 K LaMnOa (LMO) undergoes a structural 
phase transition, transforming from cubic to orthorhom¬ 
bic symmetry. Under ambient conditions the orthorhom¬ 
bic perovskite has a paramagnetic spin structure. Below 
the Neel temperature of Tn « 140 K,® LMO is an in¬ 
sulator with A-type antiferromagnetic (A-AFM) spin or¬ 
dering. In the low T orthorhombic P nma LMO shown in 
Figure [l] experimental reports of lattice parameters are 



Figure 1: (001) face of A-type antiferromagnetic (A-AFM) 
LaMnOa . Mn in purple, La in green and O in red. Arrows 
indicate direction of spin polarization on Mn ions. 


a = 5.736 A, b = 7.703 A and c = 5.540 A by neutron 
powder diffraction^. To support the A-AFM ordering 
in LaMnOa, the Mn d 4 electrons exchange anisotropi- 
cally: ferromagnetic (FM) coupling exists between Mn 
in {010} planes while AFM coupling exists between suc¬ 
cessive planes along [010]. 

The reported experimental band gaps in LMO cover a 
range of values, depending on whether the gap is deter¬ 
mined from measurements on conductivity (0.24 eV)Pj 
optical absorption (1.1 eV).® photoemission (1.7 eV)p^ 
optical conductivity (1.9 eV),® or resonant Raman spec¬ 
troscopy (2 eV).® DFT is a single-particle theory, so 
even with the exact exchange-correlation functional, it 
can only describe the fundamental (quasiparticle) band 
gap and not the optical one. We therefore consider the 
most appropriate reference value to be the 1.7 eV photoe¬ 
mission gap measured by Saitoh et al®. We note that 
recent computational work by Lee et aZ.® predicts a di¬ 
rect gap of 1.1 eV, in agreement with the optical absorp¬ 
tion gap of 1.1 eV measured by Arima et a/.®. The value 
of the optical gap is generally lower than the fundamen¬ 
tal gap due to electron-hole interactions (i.e., excitonic 
effects). Such two-particle interactions are not included 
in standard one-particle DFT, so we believe the most re¬ 
liable comparisons should be made between a benchmark 
indirect experimental photoemission gap such as the 1.7 
eV Saitoh gapt® and the indirect DFT gap. 

One of our main practical considerations here is to re- 
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Table I: LDA+? 7 e ff and GGA+[/ e ff results for the energy gap 
E Ga,p (in eV), for the A-AFM and FM phases, and the total 
energy difference A E = E A ' A¥M — E FM (in meV) per formula 
unit of LaMnOa. The crystal structure is held fixed at the 
experimental geometry. 


t/ eff (eV) 

LDA (PZ 81 ) 

GGA (PBEsol) 

77 'Gap 
^ A-AFM 


pGap 

-^A-AFM 

A?m P 

A E 

0 

0.0 

0.0 - 

-22 

0.2 

0.0 

-13 

2 

0.5 

0.0 

-5 

0.6 

0.0 

1 

4 

1.0 

0.0 

4 

1.0 

0.0 

8 

6 

1.3 

0.1 

10 

1.3 

0.1 

14 

8 

1.4 

0.2 

14 

1.4 

0.2 

17 


produce the different facets of the above experimental 
description. To do this, DFT calculations are performed 
screening through different levels of Coulombic localisa¬ 
tion. 


A. Description of LaMnC>3 using DFT+ U 

Previous work has applied the sin gle term U e g ap¬ 
proach to calculations on bulk LaMn 03 . ^ 13 l 23 l 5( ^ The fail¬ 
ure of this approach to simultaneously describe the en¬ 
ergy gap, structure and magnetism drives us to systemat¬ 
ically examine the U e g method. These initial results also 
provide context for the more sophisticated U\J method 
and analysis of its merits and behaviour below. 


1 . Experimental LaMnOa structure via DFT+f 7 e (f 

Standard LDA (PZ81) and GGA (PBEsol) with U e g = 
0 eV both successfully stabilize the low temperature ex¬ 
perimental A-AFM ordering as shown in Table [T| How¬ 
ever this is essentially where the success ends. As noted 
previously, both GGA and LDA are often unable to 
produce significant orbital splitting (beyond some as¬ 
pects due to spin exchange and structural distortion) 
and also exaggerate electron delocalization due to inexact 
exchange (or equivalently lack of self-interaction correc¬ 
tion). This inevitably results in a qualitatively incorrect 
electronic structure with a seriously underestimated band 
gap: both GGA and LDA with U e g = 0 eV yield band 
gaps that are far too small compared to experiment. 

Increasing U e g stabilizes the occupied (f ir7 > |) eigen¬ 
states and drives orbital occupations toward binary po¬ 
larization: filled states become more filled and empty 
states more empty. For example, increasing U e g from 
0 to 8 eV in GGA calculations results in the following 
change in occupancies in the Mn d manifold: 


/ 0.65 

0.22 \ 


/ 1.00 

0.10 \ 

0.73 

0.26 


0.57 

0.17 

0.93 

0.11 

-> 

0.98 

0.04 

0.93 

0.10 


0.97 

0.03 

\ 0.93 

o.o9 y 


\ 0.97 

0.04 J 


where the ordering of orbitals in the local basis is 
/ 3z 2 - r 2 \ 
x 2 — y 2 
xy 
yz 

V XZ ) 

The Hubbard limit of very large U e g typically favours 
FM coupling in LaMn 03 ,®and Table [i] confirms this. The 
primary reason is that increasing U e g kills the superex¬ 
change mechanism, which scales as ~ t 2 /U e g where t is 
the effective Mn-Mn hopping, and this mechanism under¬ 
lies the stability of A-AFM ordering in LaMnC> 3 . As a 
result, e g double-exchange is relatively strengthened and 
we find FM ordering. Critically, Table [I] shows that a 
large U e g value is required to open a satisfactory energy 
gap. Unfortunately this situation results in a trade-off 
between correct gap or correct magnetism. 


2. Relaxed LaMnOa structure via DFT+[/ e ff 

When we permit the structure of LaMnCU to fully relax 
during the calculation, we find the results in Table [Tl] We 
see that having a non-zero U e g improves the crystal ge¬ 
ometry and the electronic structure description for both 
GGA and LDA. Particular improvements are for the large 
erroneous distortion in a (insufficient orthorhombic char¬ 
acter) and the opening of the band gap. Figure[2]and Ta¬ 
ble [H] show that the band gap increases roughly linearly 
with U e g at first but then tails off at higher U e g. The 
ineffectiveness of U e g at high values is shown in Figure 
[U and can be understood in terms of the partial x 2 — y 2 
occupation shown in Equation ([H]). The partial x 2 — y 2 
occupation damps the impact of U e g on the energy eigen¬ 
value splittings since Ae x 2 _ y 2 a = U e g (| — f x 2 - y 2 a ) ~ 0 
for fx 2 -y 2 a ~ (The reason partial e g occupation oc¬ 
curs is that the d manifold is not isolated but connected 
to the rest of the system via hybridization to the O 2 p 
orbitals, or in other words due to the partial covalency 
of the Mn-0 bond.) 

Although adding U e g to GGA and LDA produces sim¬ 
ilar band gaps as per Table [TTJ the GGA+t/ e ff geometry 
is preferable. Overall a high value of U e g ~ 6 eV, cor¬ 
recting the GGA formalism, provides on balance the best 
gap/structure combination. Again, an evident failure of 
U e g is its inability to predict the correct magnetic order¬ 
ing at the U e g level required to correct the structure and 
the band gap. 
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Table II: Band gap E Gap , total energy difference A E = e a ~ afm _ E FM p er unit cell, and percent errors, with respect to 
experiment, for lattice parameters and unit cell volume of fully relaxed A-AFM bulk LaMnOs. 
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19 
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Figure 2: t/ e ff (left) and U\J (right) corrected density of states for fully relaxed A-AFM LaMnOs, as a function of energy, 
E — Ep. Black curves show the total density of states while red and yellow curves show Mn d majority spin and minority spin 
densities of states. 


3. Experimental LaMnOs structure via DFT+[/| J be verified by fully relaxing ionic positions and lattice 

parameters, which we report on next. 


Following the failure of the U e g scheme in both single¬ 
point and relaxed geometry calculations, we turn to the 
DFT+C/|J methodology. The U e g results conveniently 
suggest a reasonable starting point: since U e g = U—J, an 
U | J correction with magnitude of approximately U — J ~ 
6 eV is appropriate. Results in Table [ill] are for bulk 
LaMnOs at the experimental structure, and sample J 
from 0 to 3 eV in conjunction with U = 6 eV and 8 eV. 
Increasing J for a fixed value of U stabilizes A-AFM or¬ 
dering and enhances orbital splitting which further opens 
the band gap. Orbital splittings due to the U e g are gen¬ 
erally “isotropic” in that they are based solely on the oc¬ 
cupation. The marked improvement by the U\ J method 
emphasises the importance of explicit spatial exchange 
anisotropy in the LaMnOs Mn d manifold. The results 
in Table m are encouraging, but since distortion of the 
lattice is critical in LaMnOs P the trends observed must 


4- Relaxed LaMnOs structure via DFT+[/| J 

Table [IV] displays key data for fully relaxed bulk 
LaMnOs using the DFT+(7|J framework. Relaxed re¬ 
sults largely echo the experimental structure results 
above for A-AFM LaMnOs, with the U\J combination 
of U = 8 eV and J = 2 eV providing a good material 
description. In particular the U\J = 8|2 eV combina¬ 
tion provides agreement in terms of electr onic, magn etic 
and structural observables from experiment 44 46 47 — and 
also more computationally expensive many body GW 
approximation 52 r esults. Volume errors < 1 % improve 
on previous work, 23 44 51 and the error in energy gap 
is small at approximately ~ 5 % (< 0.1 eV errorpS 
In addition the A-AFM ordering is stabilized against 
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Table III: Results from GGA (PBEsol) + U\J for the exper¬ 
imental geometry of bulk LaMnOs. Band gaps J5 Gap are in 
eV for each magnetic state, and A E = e a - afm _ £; FM ; s the 
total energy difference per unit cell between the two magnetic 
phases. 


U\J (eV) E^£ fm (eV) (eV) A E (meV) 


6|0 

1.3 

0.1 

14 

6|1 

1.3 

0.2 

6 

6|2 

1.2 

0.4 

-10 

6|3 

0.6 

0.2 

-39 

8|0 

1.4 

0.2 

17 

8|1 

1.5 

0.4 

11 

8|2 

1.6 

0.8 

-2 

8|3 

1.0 

0.5 

-32 


FM orde ring which was previously seen as a missing 
ingredient 23 * 24 !. The improvements in LMO description 
depend intimately on the intra-orbital exchange descrip¬ 
tion - this is explored further by quantifying the action 
of the Hund’s coupling interaction on the LMO Mn d 
states. 


B. Explicit exchange anisotropy in Mn 3+ 

Strong on-site Coulomb repulsion is the central theme 
in paradigms of “Mottness” and electron localisation. 
However, the importance of Hund’s coupling (intra¬ 
orbital exchange) in materials with partial d and /-shell 
occupations has been highlighted.^ In this section, we 
attempt to understand the nature of Hund’s coupling in 
LMO, by examining the effects of the on-site exchange 
terms as defined in Appendix A. We explore why the 
U\J methodology can describe LaMnC >3 adequately, re¬ 
producing band gap and correct magnetic ground state 
simultaneously. We employ a simple model where we fo¬ 
cus only on the occupancies of the Mn 3+ d 4 manifold in 
order to isolate the effect of the exchange J parameter 

g id related physics) on the Mn d states as per Equation 
. Majority spin t 2g 3 states are generally fully occupied 
e to the strong exchange splitting between spin chan¬ 
nels, and as is well known, increasing U increases occu¬ 
pancy polarization. However the nature of anisotropic 
interactions in the Mn d shell due to J is less well un¬ 
derstood, particularly with respect to the polarization of 
the e g 4 occupation into 3 z 2 — r 2 or x 2 — y 2 (or some mix 
of the two). 

We begin with three model e g 4 occupations, n eg = 
0, ±1, in an attempt to pinpoint what J really does and 
understand the nature of Hund’s coupling in different 
limits. As a reminder, 7r e s is the e g occupancy polariza¬ 
tion as defined in Equation 0. After examining these 
model systems, we will consider the effect of J in the 
actual calculations where we use the calculated ab ini¬ 


tio occupations together with our analytical rewriting of 
the U\J energy function and eigenvalue corrections. As 
explained above, the e g and f 2g group terms discussed 
correspond to the local octahedral basis (ie., post rota¬ 
tion as per Appendix B). 


1. Anisotropic exchange for model orbital occupations 

To illustrate the anisotropic effects of the J terms in 
the U\ J schema, we begin with a set of model occupancies 
where we fix the formal occupation of Mn 3+ (d 4 ) but vary 
the orbital polarization. 

A 7r e s = +1 model polarization corresponds to a single 
hole on the majority spin 3z 2 — r 2 site (i.e., / x 2 _. y 2 0 . = 
/t 2 gcr = H f 3 z 2 -r 2 a = 0 and fa = 0). Based on Equation 
Q, the added effect of the exchange J terms is to create 
additional energy splittings (beyond simple occupancy 
polarization proportional to U e g) given by 


7r e * = +1 : 


(Ae CT |Aeg) = J ■ 


/ 0.00 

-1.14 \ 

0.52 

0.63 

0.52 

0.63 

—0.52 

-0.06 

\ -0.52 

-0.06 / 


The opposite polarity, n e « = — 1, corresponds to a single 
hole on the majority spin x 2 — y 2 site (that is, / 3 2 2 _ r 2 CT = 
ft 2 go- = 1, f x 2 — y 2 a = 0 and fa = 0). This results in the 
following exchange energy splittings 


7T e e = -1 


(Ae CT |Aeg) = J ■ 


/ 0.52 

0.63 \ 

0.00 

-1.14 

-0.86 

-0.29 

0.17 

0.40 

\ 0.17 

0.40 / 


Removing the polarization, 7r eg = 0, the single hole is 
equally spread over the two majority spin e g sites (/ eg(T = 
0.5 and /t 2g0 - = 1 and fa = 0). This leads to the splittings 


7r e s = 0 : 


(Ae CT |Ae ff ) = J • 


/ 0.26 

-0.26 \ 

0.26 

-0.26 

-0.17 

0.17 

-0.17 

0.17 

\ -0.17 

0.17 / 


( 8 ) 


To visualize these results, we display a schematic showing 
these splittings in these three cases of 7r 6g = 0, ±1 in 
Figure [3] where the corrections due to both U and J are 
shown. 

These model results together with the Figure [3] clearly 
point out that the effect of the J terms is explicitly 
anisotropic and its anisotropy and precise value depends 
on the orbital polarization (which may have been present 
due the action of the U e g term). The anisotropy ex¬ 
ists across both the magnetic quantum number and spin 
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Table IV: Fully relaxed LaMnOa results based on GGA (PBEsol) + U\J. Band gaps if Gap , lattice parameter errors, and total 
energy differences between the A-AFM and FM magnetic phases A E = £ A ~ AFM _ £ FM per formula unit are listed. 


U\J (eV) 



A-AFM 



A E (meV) 

£ Gap (eV) A a (%) A b (%) A c (%) A 

VfaL (%) 

6|0 

1.1 

0.4 

-0.1 

-0.4 

0.0 

19 

6|0.5 

1.2 

0.6 

-0.4 

-0.4 

-0.2 

14 

6|1 

1.3 

0.9 

-0.7 

-0.5 

-0.3 

7 

6| 1.5 

1.4 

1.2 

-0.9 

-0.6 

-0.3 

-1 

6|2 

1.4 

0.8 

-1.2 

0.1 

-0.4 

-10 

6|2.5 

1.3 

1.8 

-1.3 

-0.8 

-0.2 

-21 

7|0 

1.2 

0.5 

0.1 

-0.3 

0.3 

23 

7|0.5 

1.3 

0.7 

-0.2 

-0.4 

0.2 

17 

7|1 

1.4 

0.9 

-0.5 

-0.5 

0.0 

11 

7| 1.5 

1.5 

1.2 

-0.7 

-0.5 

-0.1 

3 

7|2 

1.6 

1.6 

-1.0 

-0.7 

0.0 

-6 

7|2.5 

1.5 

1.9 

-1.2 

-0.7 

-0.1 

-17 

8|0 

1.1 

0.5 

0.4 

-0.3 

0.6 

27 

8|0.5 

1.2 

0.6 

0.1 

-0.3 

0.4 

21 

8|1 

1.4 

1.0 

-0.2 

-0.5 

0.3 

15 

8| 1.5 

1.6 

1.2 

-0.5 

-0.5 

0.1 

7 

8|2 

1.8 

1.5 

-0.8 

-0.6 

0.1 

-2 

8|2.5 

1.7 

1.7 

-1.0 

-0.6 

0.1 

-14 


channels (cr and a). We now discuss these three cases in 
more detail. 

When we have full e g orbital polarization, i.e. 7r e ® = 
±1, each polarity produces a unique splitting pattern 
where the magnitude of anisotropy depends on the sign 
of 7 r eg . This is despite the fact that 3 z 2 — r 2 and x 2 — y 2 
states both have e g symmetry: as we can see that occupy¬ 
ing each one (separately) splits the f 2g quite differently. 
This difference is due to fact that the x 2 — y 2 state is sym¬ 
metry related to the f 2g states (it is the xy state rotated 
by 7 t/ 4 about the z axis). For example, when x 2 — y 2 
is fully occupied, the splittings for x 2 — y 2 and xy are 
identical but differ from the other orbitals, but the same 
is not true when 3 z 2 — r 2 is filled instead. Interestingly 
J2icr Aeiafia = 0 when 7r eE = ±1: this indicates that 
neither polarization is energetically preferred by intra¬ 
orbital J terms. 

With zero e g orbital polarization, i.e. ir £s = 0, we 
find that this degeneracy inhibits anisotropy from the 
J terms: the splitting within each f 2g and e g manifold 
is isotropic for both spin channels. The action of the 
J terms in this situation is to shift the energies of this 
manifold en masse. Here, Ylicr A e i(J f i(J = — 0.26 J when 
7 r e ® = 0 compared to zero for 7r e « = ±1. Hence, the 
anisotropic exchange terms in isolation actually favour 
degenerate occupancy. This result appears to be counter¬ 
intuitive given the importance of J to anisotropy. The 
resolution is that we have a much larger and dominant di¬ 


rect Coulomb term U that produces orbital polarization 
in the first place; the weaker J terms then further enlarge 
the polarization and make the system more anisotropic. 
Table |V| shows this behaviour numerically. 

In brief, we see that J acting alone energetically 
favours degeneracy. However, with a strong U term al¬ 
ready creating orbital polarization, the J terms provide 
the enlarged anisotropic splitting that one finds in the 
final results of the calculation. 


2. Anisotropic exchange for ab initio orbital occupations 

For the ab initio orbital occupations we use the 
DFT+[/| J = 8|2 eV calculation results, which yield an 
occupation-polarized e g manifold as shown in Table [V] 
The e g polarity is found to be orbitally ordered across 
the LaMnC >3 unit cell as shown in Figure [4] We now 
examine this situation in more detail. 

The Mn d occupancies from the U\J = 8|2 eV method 
with relaxed geometry are 


/ 0.99 

0.10 \ 

0.41 

0.33 

0.98 

0.06 

0.96 

0.04 

V 0.97 

0.05 / 
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Figure 3: The occupation of states is represented in the model Mn 3+ d 4 manifold (majority spin only). Orbital degeneracy is 
broken by octahedral crystal field (CF), Coulomb repulsion U e a (U in the figure) and exchange J following Equation Q. Each 
vertical bar represents one unit of electron occupation. 7r e9<T is defined in Equation 0 for which three limits are examined: 


= 0 {fx 2 -y 2 a = /3z 2 -r 2 a = 0.5), 


= + 1 {fx 2 -y 2 a = f3z 2 —r 2 a + 1 = 1), and = -1 (f x 2 _ y 2 a + 1 = f 3z 2 _ r 2 a = 1). 


Local octahedron basis ^ Mn e g orbital in local basis 





Short basal , 

Mn-0 / 
x 


y Apical 
I Mn-0 


0.99 3z 2 -r 2 


0.41 x 2 -y 2 1.97 z 2 -0.58 x 2 -1.4 


Long basal 
Mn-0 
z^' 
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( 010 ) 


( 001 ) 


Figure 4: Orbitals in the LaMnOs unit cell from a C7| J = 812 eV calculation, a) MnOe octahedron with Jahn-Teller distorted 
plane and local basis vectors labelled, b) Visualization of the occupation of the 3z 2 — r 2 and x 2 — y 2 states in the rotated basis 
of the density matrix as well as their superposition for the total local e g occupancy (plotting the occupation times the orbital 
expressed in spherical harmonics), c) The ordering of the occupied e g shell (1.97z 2 — 0.58a: 2 — 1.4 y 2 ) in the (010) Jahn-Teller 
distorted FM coupled plane, d) The ordering of the occupied e g shell in the (001) plane with AFM coupling along b. Note 
x, y, z is the local octahedron basis, and a, b, c lattice vectors correspond to x', y', z’ global (pre-rotation) calculation basis. 
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These occupancies correspond to 7T egCr = —0.41 (from 
/ 3 . 2 _ r 2 CT = 0.99 and f x 2 - y 2 a = 0.41). The Mn d shell ob¬ 
viously has more electrons than the model system above 
which was based on formal occupancies for Mn 3+ . Again, 
we note that this increase is due to hybridization of the 
Mn d orbitals with the neighboring O 2 p orbitals which 
admixes some Mn d into the low-energy occupied va¬ 
lence states and increases the electron count. Put differ¬ 
ently, itineracy due to the kinetic energy minimization 
competes with Hubbard-esque Coulomb repulsion and we 
reach a balance. Numerically, for U\J = 8|2 eV, the oxi¬ 
dation state based on the above Mn d occupations can be 
calculated to be 2 . 12 + (an alternative or complementary 
Bader charge picture yields an oxidation state of 1.68+, 
still less than the formal 3+). 

The U | J = 8|2 eV occupancies of Equation ([9]) result 
in energy splitting beyond splitting from U e s alone: 


(Ae CT |Ae ff ) = J ■ 


/ 0.15 

0.20 \ 

0.30 

-0.65 

—0.41 

-0.06 

-0.03 

0.18 

\ - 0.01 

0.20 J 


( 10 ) 


The e g occupancy polarization of the XJ\J = 8|2 eV cal¬ 
culation is considerably weaker than the previous model 
cases. Nevertheless, it is large enough to drive significant 
anisotropic exchange splittings in Equation (10 1 . For ex¬ 
ample, the splittings are anisotropic within the majority 
spin f 2 g manifold: the xy state is pushed down by ap¬ 
proximately 0.4J compared to the other two f 2 g states. 
Within the e g manifold, the fully occupied 3 z 2 — r 2 state 
is pushed up by 0.15J while the partially occupied x 2 —y 2 
state is pushed up considerably more by 0.30J. 

The direct Coulomb interaction U obviously increases 
7 r 6g , as expected from the standard instability condition 
for orbital polarization 2 ^!!, 


U eS x D”(E ¥ ) » 1, 

where D cr (Ep) is the density of states in the a spin chan¬ 
nel at the Fermi level. The origin of the monotonic rela¬ 
tion between 7 r 6g and J, shown in Table [Vj is less obvious 
as J is naively expected to drive the electronic structure 
away from orbital polarization as U e g = U — J. However 
7 r 6g does increase with J, for the above-noted reason that 
J alone may favour orbital degeneracy but J is strongly 
anisotropic when in conjunction with a large U value, 
resulting in the unequal upward “push” of the two e g or¬ 
bitals with increasing J. That J and 7 r eg are so strongly 
coupled in this material is interesting, as band gap, Jahn- 
Teller distortions, and inter-site magnetic couplings all 
depend on 7r eg . 

As first examined by Kugel and KhomskiiP^ e g 1 oc¬ 
cupation polarization (i.e., an electron-electron Jahn- 
Teller degeneracy breaking) enhances virtual superex¬ 
change interactions, relative to kinetic exchange inter¬ 
actions such as FM double-exchange. This competition 


Table V: Orbital occupation polarization, 7r egCT , for fully re¬ 
laxed LaMnOa within the U e b and U\J approaches. Majority 
spin are a and minority spin are a. 


Correction (eV) 

Like-spin 
polarization, n e9<T 

Opposite-spin 
polarization, 7r es<T 

U e B = 0 

0.06 

0.08 

aj 

II 

00 

-0.27 

0.25 

a 

II 

00 

H 

-0.33 

0.41 

CM 

00 

II 

A 

-0.41 

0.54 

CO 

00 

II 

"-s 

A 

-0.52 

0.65 


between superexchange and double-exchange is observed 
in the LMO magnetic ground state, which varies accord¬ 
ing to the magnitude of 7 r eg (e g occupancy polarization). 
7 r 6g increases with J, which explains the flip in long range 
magnetic ordering of the ground state from FM to A- 
AFM as the intra-orbital parameter J is increased. 

At J = 2 eV the value of 7 r eg in Table |V| is large enough 
to stabilize the correct A-AFM ordering (as per Table 
IV). The U\J = 8|2 A-AFM ground state in Figure [l] 


corresponds to a 0.99(3z 2 — r 2 ) + 0.41(;r 2 — y 2 ) e g occu¬ 
pation density in the local octahedral basis. The orbital 
ordering pattern across the unit cell is shown in Figure 
[4j and can be rationalized in t erms of the Goodenough- 
Kanamori superexchange rules 1SS1SZ1 

The 0.99(3z 2 — r 2 )+0Al(x 2 —y 2 ) e g occupation density 
can be rewritten as 1.97 z 2 — 0.58a: 2 — lAy 2 . This expres¬ 
sion shows the anisotropy in the e g state, in particular 
between the z and x directions in the octahedron: the 
z 2 contribution is much larger than a: 2 , as per Figure [4] 
b. Each octahedral frame in the ac plane is related to its 
neighbour by a 7r/2 rotation about the b lattice vector, 
so z 2 /x 2 anisotropy forms a checker board pattern of e g 
partial occupation in the ac plane. Note this corresponds 
to the long/short Jahn-Teller Mn-0 pattern in the ac 
plane, as per Figure[4]c). According to the Goodenough- 
Kanamori rules, superexchange in the ac plane is deter¬ 
mined by z 2 / x 2 anisotropy in the e g partial occupation, 
and results in the FM coupling in the ac plane. 

The y 2 component of the e g partial occupation forms 
occupied stripes pointed along local octahedra y axes, 
following the b lattice direction (with a small tilt) as 
in Figure [4] d). The continuous stripes of y 2 character 
along the b lattice direction correspond to the ’non-Jahn- 
Teller’ Mn-0 bonds in this direction. The Goodenough- 
Kanamori rules determine that the continuous stripe of 
y 2 character from the e g partial occupation corresponds 
to AFM superexchange. The AFM coupling is along the 
b lattice parameter direction, between the FM coupled ac 
planes. Together the in-plane FM and inter-plane AFM 
produce the A-AFM ground state of LaMnOa, so our 
U\J = 8|2 eV calculation results are in-line with experi¬ 
ment as well. 

If instead the Hund’s coupling was weaker (i.e., smaller 
J), then 7 r 6g would also be smaller. This alters the char- 
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Figure 5: LaMn 03 densities of states (DOS) for U\J = 8|2 
eV and U\J = 813 eV calculations. Majority spin corresponds 
to positive DOS and minority to negative DOS. 


acter of the occupied states in the e g shell, so that orbital- 
ordering mediated A-AFM superexchange is reduced rel¬ 
ative to other effects such as FM double-exchange. This 
explains why stabilization of A-AFM magnetic ordering 
(see Table IV) is only possible when intra-orbital ex¬ 
change is large enough. Too small of an intra-orbital 
exchange interaction is the origin of the incor rect FM 
ground state found in prior examination^ 2 23 EH of LMO 
using standard DFT. 

The improvements in the LMO description through ap¬ 
plying exchange corrections reinforce hints by Sawada et 
fi/.pS Solovyev et a/P and Hashimoto et oZ.J^that the 
correct orbital and magnetic ordering in LMO requires an 
anisotropic intra-orbital exchange correction to the DFT 
ground state. In what follows, we discuss further details 
of the electronic and crystal structure. 



Figure 6: Local magnetic moment within the Bader volume 
for Mn cations and O anions (within the (010) basal plane), 
from DFT+t/eff and DFT+t/| J methods. The notation U§\J 
and Us \ J indicates U is fixed to 6 eV and 8 eV, respectively, 
while J is varied. The experimental reference local magnetic 
moment is 3.7 


high sensitivity of the electronic structure of LMO to 
perturbations in part underlies its complex phase dia¬ 
gram. This is illustrated by comparing the XJ\J = 8|2 
and U\J = 8|3 eV DOS in Figure [5j and examining the 
magnetic state of DOS near the edges of the valence 
band maximum (VBM) and conduction band minimum 
(CBM). For U\J = 8|2 eV, Hund’s rules are obeyed as 
both VBM and CBM have the same spin state whereas 
Hund’s rules are broken for U\ J = 8|3 eV. We find that 
the cross-over occurs at 2.4 eV. LaMnOa is fragile in 
terms of exchange: above J ss 2.4 eV we have the break¬ 
down of Hund’s rules while below J ss 1.8 eV incorrectly 
stabilizes the FM rather than A-AFM ground state. 


C. Electronic and crystal structure details 


2. Magnetic coupling constants 


1. Orbital order 


It was previously shown that applying Coulomb cor¬ 
rections, such as with XJ\J = 8|2 eV, corrected the 
LMO electronic, magnetic and lattice structure. Fur¬ 
ther electronic structure details are shown for the LMO 
DOS in Figure [5| at the U\J = 8|2 eV level of correc¬ 
tion. In Figure]!)] the position of each band in the Mn 
DOS agree quantitatively with the optical conductiv¬ 
ity measurements of Jung et al.. i7 Further experimen¬ 
tal agreement comes from our U\J = 8|2 eV calculated 
local magnetic moment, which at 3.7 /tb agrees with El- 
eman’s measurement EQ The U and J dependence of the 
local magnetic moments are shown in Figure [bj The 


The magnetic coupling constants in LMO have b een 
extracted by Munoz et alW^ amongst otherj^ 4 * 49 * 58 ( by 
considering an Ising model (with S = 2 spin moment per 
Mn ion) for different spin-ordered solutions. The intra¬ 
plane ( ac) J\ and inter-plane (b) Ji coupling constants 
for the 20 atom LMO unit cell are 


= sW 
* = 


G-AFM _ ^A-AFIvn 


A-AFM _ ^FMj 


( 11 ) 


The initial A-AFM/FM stability results in Table IV 


hint that the coupling constants will depend strongly 
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on the Hund’s exchange parameter. In the context of 
previous works, Ji are well known to be highly sensi¬ 
tive, for example to variation in Mn-O-Mn angle through 
superexchange interactional, and the Mn ionic charge 
population^. 

On the trend of magnetic stability in U and J, the su¬ 
perexchange interaction which stabilizes AFM ordering is 
expected to scale inversely with effective on-site Coulomb 
interaction, i.e. ~ t 2 /U e s where t is effective inter-site 
hopping. Considering first J 2 (~ E A ~ AFM — E FM ) in Fig¬ 
ure [ 7 ] the stability of AFM coupling along b decreases 
both with increasing U or decreasing J as expected since 
Ueg = U — J. However, the dependence of the J 2 cou¬ 
pling on U and J is not equivalent: the variation in J 2 
is some five-fold more sensitive to changes in J than U, 
i.e., <9J 2 /<9J ~ —5 dJ^/dU . The AFM coupling in the ac 
plane, measured by ,J\ , is even more sensitive to the intra¬ 
orbital Hund’s parameter, with dJi/dJ ~ — 10dJi/dU. 
The origin of coupling constant sensitivity to J, is the 
strongly anisotropic effect of J on the LMO Mn d shell 
states, with variation in J increasing 7r e » in Table |V| 
above and beyond that accesible with U alone. 

Neutron scattering experiments have determined cou¬ 
pling constant values of J? xp = 0.83 meV and J ^ xp = 
—0.58 meVP3 In Figure [7j reasonable values for J\ are 
produced with J ss 1.75 eV, and for J 2 with J 2 « 2 eV. 
The discrepancy in J value for each Ji is perhaps unsur¬ 
prising given the extreme sensitivity of A-AFM, G-AFM 
and FM phases to the intra-orbital Hund’s interaction. 
Overall the U\J = 812 eV combination previously sug¬ 
gested remains a good compromise at the level of half 
integer eV screening intervals considered here. Although 
higher resolution screening in J is beyond the scope of 
this work, if DFT+(7| J calculations are required for ther¬ 
modynamics applications, results indicate a small modifi¬ 
cation of J by a few percent may be advantageous to tune 
the magnetic transition temperatures precisely, while the 
magnetic couplings are relatively insensitive to the direct 
U term. 

Due to the sensitivity of the magnetic couplings to the 
Coulombic J correction, agreement with experiment is 
challenging. At the U = 8 eV required to open the band 
gap, and screening in J at half integer intervals shown 
in Figure [7| U\J = 8|2 eV remains the best compro¬ 
mise. For XJ\J = 8|2 eV the inter-plane coupling at 
J 2 = —0.30 meV has the correct sign but in magni¬ 
tude falls short of = —0.58 meV®. More problem¬ 
atic is the intra-plane coupling, which overestimates the 
tendency for electrons to couple antiferromagnetically in 
the ac plane, excessively stabilizing G-AFM ordering at 
J 2 = —0.19 meV compared to J| xp = +0.83 meV. 

In Figure [ 7 ] the colored areas show the J values that 
correspond to coupling constants between zero and J® xp , 
i.e. the correct sign for each Ji. The overlap in col¬ 
ored areas identifies the narrow range of intra-orbital 
exchange values, 1.88 < J < 1.95 eV, that gives the 
correct signs for both Ji together, with A'(A-AFM) < 
A'(FM) < E(G- AFM) in agreement with experiment®. 



Figure 7: LaMnOa magnetic coupling constants Ji and J 2 
versus U\J schema Hund’s exchange parameter J, for U = 6 
eV (white squares), U = 7 eV (white circles), (7 = 8 eV (black 
squares). Ji and J 2 are defined in Equation (111. The red- 
blue overlap (centered at J = 1.9 eV) suggests a J exchange 
value for the U\J scheme that provides the correct sign for 
both coupling constants - see main text for discussion. 


Based on the refinement in J value, we have performed 
GGA (PBEsol) calculations with U\J = 811.9 eV. The 
error in calculation results with respect to experimen¬ 
tal valued^] 47 51 l 6 °l tSH j s summarized in Table lvil The 
U\J = 8| 1.9 eV calculations produce good experimen¬ 
tal agreement overall, with magnetic coupling constants 
with signs that agree with experiment, Ji = +0.2 and 
J 2 = —0.1 meV, a band gap value only a couple of per¬ 
cent above the experimental 1.7 eV value, lattice param¬ 
eter errors between +1.5 % and —0.8 % which largely 
cancel to give a volume error with respect to experiment 
of +0.1 %™ 


3. Jahn-Teller distortion 

We end our analysis with conclusions on the nature of 
Jahn-Teller distortion in LMO and on the origin of the 
LMO insulating state. Jahn-Teller distortion in LaMnOa 
can be characterized in terms of two normal modes of the 
type introduced by van Vlecli® and by Kanamori®. The 
normal modes are shown in Figure [8] along with the crys¬ 
tal unit cell and the local octahedral basis. The modes 
are calculated as 

Q° rtho = [Y 2 - Y 5 - Xx + X 4 ] 

Q Tetra = -4 [2Z 3 - 2Z 6 - Y 2 + Ys -X , + X 4 ] . 

v 6 

Each variable represents an octahedral bond length, with 
subscripts indexing oxygen octahedral cage sites for a 
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Table VI: LaMnOa electronic, magnetic and structural properties obtained from a U\ J = 8| 1.9 eV calculation, with comparison 
to experimental counter-parts.^ 47 ! 51 1 60 ^^ The J = 1.9 eV value is based on a refinement of the Hund’s exchange parameter to 
secure the correct sign for both magnetic coupling constants, J i and J 2 , which are exceptionally sensitive to on-site exchange 
- see Figure [7] 

,,, j , ... Electronic gap Magnetic Structural 

1 (C ’ E Gap (eVj 13 CharacteJ^ J u J 2 (meV^ M Q° rtho , Q Xetra (a.n.jSl a> c (Ag< V (A 3 j^> 

8| 1.9 1.75 e'g-y -+ e 2 g ^ +0.2,-0.1 3.76 0.145,0.856 5.823,7.642,5.508 245 

Exp. 1.7 ej f ef n +0.83,-0.58 3.7 0.14,0.78 5.736,7.703,5.540 245 


given manganese centre, i(i = 1, ...6). In the local basis 
in this work, which differs from other clioiceiP^M^l, Z., = 
z o _ jjMn are j Qn g ]y[ n _Q bonds and X; = xf — 
x Mn short Mn-0 bonds, with both in the FM coupled ac 
plane. Y; = yp — y Mn are along the inter-plane AFM 
coupled b lattice direction. 

We begin commenting that the formation of a band 
gap in LaMnC >3 is not solely electron-electron (e-e) 
or electron-lattice (e-l) in character. Rather, it is a 
joint function of the lattice relaxation and development 
of Jahn-Teller distortions as well as the strong on-site 
Coulomb interaction. This is illustrated explicitly in Fig¬ 
ure [8| As mentioned above, two logically distinct routes 
to breaking symmetry exist in order to produce a gap: 
(i) a purely electronic effect via electron-electron interac¬ 
tions and the formation of a sizable orbital polarization 
7 r e ® that breaks symmetry (also called e-e Jahn-Teller 
distortion]^, or (ii) electron-lattice (e-l) Jahn-Teller dis¬ 
tortions where certain local octahedral phonon modes be¬ 
come soft, the Mn-0 bond lengths become unequal, and 
this creates crystal field symmetry breaking. These two 
mechanisms are in fact mutually enhancing, and which 
one causes which in LaMnOa is an open question that has 
been debated in the works of KhomskfP^, Yin et al 
and Loa et a m. 

In some materials, one mechanism can clearly domi¬ 
nate over the other. For example, in KC 11 F 3 , to which 
LaMnOa is superficially similar as both are perovskites 
with partial e g occupation, the symmetry lowering is 
truly driven by electronic interactions alonej^l and thus 
KCUF 3 is said to exhibit e-e Jahn-Teller distortion. Fig¬ 
ure [ 8 ] shows that the nature of Jahn-Teller is different in 
LaMnOs. 

Firstly, with the e-l distortion frozen out, one can gen¬ 
erate symmetry breaking and a gap for a Coulomb inter¬ 
action strength (U — J) above a critical value ~ 2 eV, 
so that in principle, the lattice distortion is not neces¬ 
sary to create a gap. However, in practice, the gap and 
orbital splitting remain small without lattice Jahn-Teller 
distortions. Secondly, with U — J set to zero, the DFT 
calculations do produce weakly active e-l Jahn-Teller dis¬ 
tortions of Q° rtho = -0.02 a.u. and Q Tetra = 0.14 a.u., 
but the gap remains essentially zero. The addition of 
Coulomb repulsion via U e g greatly enhances the e-l dis¬ 
tortion of each mode to approximately Q° rtho ss 0.12 
a.u. and Q Iotra « 0.62 a.u.. However even with U e g 
applied Q° rtho and Q Tetra remain still short of experi¬ 


ment by some 13 % and 20 % respectively. As per Table 
\V\ and Figure [8] one can only go so far with U e g: the 
orbital polarization 7r e » is too weak and the Jahn-Teller 
e-l distortion remains largely unchanged with increasing 
U eS . 

The only way to bridge the deficit is through the use 
of a dedicated exchange term via the U\J approach. 
As shown in Table [Vj J increases 7r e « and anisotropy 
throughout the d manifold significantly. By increasing J 
in the U\ J scheme, the LMO Q° rtho and Q Tetra modes 
can be tuned to agree with experiment by accessing ad¬ 
ditional e-e Jahn-Teller activity otherwise unavailable. 


IV. CONCLUSION 

An isotropic Hubbard correction, such as the U e g 
methodology, is unable to simultaneously reproduce the 
band gap, experimental level of Jahn-Teller distortion 
and magnetic ordering of bulk LaMnOs. At small U e s, 
A-AFM magnetic ordering is correctly stabilized but the 
gap and structural distortions are underestimated. With 
increasing U e g values, the gap and crystal structure are 
reproduced but FM ordering is incorrectly stabilized. 
The U | J approach, with its explicit exchange dependence 
on orbital symmetry, provides a better picture of elec¬ 
tronic, magnetic and structural properties of LaMnOs. 
The origin of the U\ J success is the Hund’s coupling ac¬ 
counted for by the spatial/orbital dependence of the ded¬ 
icated exchange terms that depend on J. These terms 
selectively polarize orbital occupation through highly 
anisotropic energy splitting in the Mn d manifold. Only 
the addition of J terms, rather than crystal field or direct 
Coulomb [/, can provide appropriate and large enough 
anisotropic splitting within the f 2 g and e g manifolds. 
Orbital order due to the short range J makes possible 
the combination of long-range FM exchange in the (010) 
plane, and AFM exchange between {010} planes, to sta¬ 
bilize the A-AFM ordered ground state. Soft phonon 
modes ( e-l Jahn-Teller) and electronic occupation polar¬ 
ization (e-e Jahn-Teller) contribute jointly to the insulat¬ 
ing state, with the latter predominant. The experimental 
Jahn-Teller distortion magnitude can only be achieved 
by adding the anisotropy J provides on top of the direct 
Coulomb occupancy polarization. The best description 
of LaMnC >3 is achieved within the PBEsol+ U framework 
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Figure 8: Left top panel: Band gap of LaMnC >3 versus U — J within the U e ff approach for fully relaxed structures where both 
electron-lattice and electron-electron interactions are active (dashed line, e-l and e-e) and for structures with the Jahn-Teller 
distortion frozen out so only electron-electron interactions are active (solid line, e-e only). Left bottom and middle panels: 
Jahn-Teller normal modes versus U — J within the U e a approach (white circles), and within the U\J approach (black circles) 
for U fixed to 8 eV and 6 eV with J varied. Right panels: Orbitally ordered and strongly Jahn-Teller active ac plane, the local 
basis convention, and Q° rtho and Q Totra moc j es . 


when U = 8 eV and J = 1.9 eV. 
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Appendix A: DFT-|- U expressions 

We begin with the U\ J rotationally invariant DFT+t/ 
total energy expression 22 written for a single atomic site 
(since the corrections are linear sums over atomic sites), 
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Ed FT+u \j = -Edft + Ejj — -Edc • 

i?DFT is the total DFT energy using some flavor of ex¬ 
change and correlation, the Coulombic +U correction en¬ 
ergy is 


Ejj = - Y^ (mam” a'\V\m'am'" a')x 

a,a' ,m z 

/ a <j' _ a <7 c \ 

\rm , mrm , "m" rm"'mrm'm"^ (7<J ') 

and the double-counting correction Eg c is 

Edc = Y {U 2 J) n °( N ° - !) + . 

(7 

In the above expressions, V(r, r') = 1/1r — r'| is the 
bare Coulomb interaction, a labels spin where a is the 
opposite spin to a, m labels angular momentum states of 
the atomic shell under consideration (d orbitals for Mn 
in this paper), is the single-particle density matrix, 

N a = trace(p a ) = J2 m ,m' Pmm'^rnm' is the number of 
electrons on the site of spin u, and U and J are the 
direct and exchange Coulomb parameters. 

The matrix elements of V are de¬ 
fined by (toctto / V'|D|to , <t77i , "c7') = 

f dr f dr' </4 CT (r)^ m v(r) 4>* m „ CT ,(r , )d> m '" (T '(r'). 

The matrix elements of V are further decomposed for an 
atomic shell with angular momentum l via 


(mam”<j'\V\m!am'"a') = 

21 

Y ° k ( lm > lm')c k (lmlm")F k 

k—0 

where c k and F k are standard atomic Slater angular 
integrals and radial integrals. For d shells, U = F°, 
J = (F 2 + F 4 )/14 and F 4 /F 2 = 0.625 are the canonical 
choices 22 } Thus only two parameters are needed to spec¬ 
ify the radial integrals: F° = U, F 2 = (112/13 )J and 
F 4 = (70/13) J. 

To make progress with expressions for Ejj and Eg c 
which are given in terms of Pm m ' and N a , we need rewrite 
these expressions in terms of the occupancy eigenvalues 
of the density matrix, f la . Denoting the eigenvectors of 
Pmm' as V m,i S0 that 


o a , = v * 7 f (v CT , Y 

rmm / j m,i jig \ m ,i) 


we may insert this expansion into the expression for Ejj. 
After some algebraic manipulations, using the fact that 
c°(lm,lm') = 6 m m' for the k = 0 term and the unitarity 


of the eigenvector V CT matrices, we find 


Eu = ^(N 2 -Yfl' S j + 


where N = N a is the total electron count on the site, 
and the Coulombic C' 7 ' 7 and exchange X CT correction 
matrices are given by 


C.f = I> 


21 

“ x 

k—2 mm'm"m n 


&m—m',m"' — m" ^ 


(V°)l m c k (lm, Im'yV^ n , i (V a ' ) c k (lm '", Zm")V^. 
and 


21 


y = Y pk Y 




k—2 


mm'm" m" 


(V^Lc fc (im,im')V", 7 . (V^L mC k (lm"' 


The Coulomb correction C <7<T matrices have zero average 
over all entries, a fact easily shown by using some basic 
properties of the Slater angular integrals. The same can 
be done for the exchange correction matrices by separat¬ 
ing out a constant term 


X£ = AX/ 7 + J( 1 - S tj ). 

Substituting this into the previous Ejj expression and 
subtracting the double-counting term E & c to cancel com¬ 
mon terms then yields the total energy 


S’dft+c/ij = -Edft H-x — ^(/i< t - fD+ 




which is in the form of the DFT+t/ e ff (Dudarev) energj^S 
plus a correction involving the C' 7 ' 7 and AX' 7 matri¬ 
ces and the occupancies. Therefore, the U\J scheme can 
be viewed as a correction to the U e g method which in¬ 
cludes additional Coulombic and exchange terms stem¬ 
ming from exchange integrals between different orbitals: 
this is because both C <7<7 and AX a are proportional to J 
and thus the orbital shape dependence of the Coulombic 
interactions on the site, something neglected by the U e g 
scheme. 

The correction to the eigenvalue follows from the occu¬ 
pancy derivative of the added terms to the DFT energy 


d(E v - E dc ) 


(U-J) + 
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In what follows, it is more convenient to work with vec¬ 
tors and matrices. Thus if we collect all occupancies f ta 
into a column vector f a , then the above eigenvalue cor¬ 
rection can be more compactly written as 

V/„ (Ejj - E dc ) = (U - J) Q - + 

J [A°f a + B CT /a] 

where we have peeled off the constant J and also in¬ 
dicated same spin and opposite spin occupancy depen¬ 
dences via the unitless matrices 

A a = _ AX CT )/J 

and 

B CT = C a9 /J. 

We now proceed to actual example cases to compute 
numerical values for the A CT and B CT matrices. The sim¬ 
plest assumption is to take the spherical harmonic states 
Yi m as the eigenbasis of the density matrix p a . This 
means V CT = I and one can directly compute the ma¬ 
trices using numerical values for Slater angular integrals. 
The results are 



( 

Y 22 

^21 

^20 

Y 2i - 1 

*2,-2 \ 



y 22 

0 

-0.52 

-0.52 

0.17 

0.86 1 


A CT = 

y 21 

—0.52 

0 

0.52 

-0.17 

0.17 



T 20 

-0.52 

0.52 

0 

0.52 

-0.52 



l F2 ’- 1 

0.17 

-0.17 

0.52 

0 

-0.52 | 



V *2,-2 

0.86 

0.17 

-0.52 

-0.52 

0 ) 



and 



However, this basis is not the most relevant for solid state 
systems such as perovskite oxides. For high symmetry 
situations, the eigenbasis of the density matrix will be 
given by t 2g (a :y, yz, xz) and e g (3 z 2 — r 2 ,x 2 — y 2 ) states. 
The conversion matrix is 


/0 l/s/2 i/V 2 0 0 \ 

0 0 0 -i/V 2 l/y/2 

V CT = 1 0 0 0 0 

0 0 0 i/V 2 1/V2 

\0 l/v/2 -i/y/2 0 0 / 


if we choose the order (3 z 2 —r 2 ,x 2 — y 2 , xy, yz, xz). The 
transformed matrices are now in the more useful basis 
with entries 


and 



These matrices directly tell us how the U\J scheme 
corrects the energy eigenvalues beyond the U e g energy 
shift. For example, the diagonal entries of B CT indicate 
that occupying any orbital pushes up the energy of the 


opposite spin orbitals by 1.14J. 

As another example, if we have an ion such as Mn 4+ 
with a full up spin and empty down spin t 2 g shell, 
so that = (0,0,1,1,1) and /j. = 0, then for the 
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up spin orbitals the energy correction beyond U e s is 
(0.52,0.52,—0.34,—0.34,—0.34) J which destabilizes the 
same spin e g and stabilizes the same spin f 2g while for 
spin down orbitals the situation is exactly reversed with 
energy correction (—0.52,—0.52, 0.34,0.34,0.34) J. A fi¬ 
nal example is a full f 2g 6 shell such as Co 3+ which gives 
zero correction to the U e s scheme. The above two matri¬ 
ces form the basis for various analyses in the main text. 


positions: = r (Oi) — r(Mn). These vectors are then 

normalized and define the local axes for the Mn. We 
create a 3 x 3 rotation matrix R connecting the global 
x' , y' , z' and local x, y , z systems 

( X '\ ( X \ ( R 11 R 12 R 13 \ / X \ 

y' I = R y = R 2 i R 22 r 23 y 
\z' ) \ Z ) \ R 31 R 32 R 3 3 J \ Z J 


Appendix B: Density matrix rotation to local axis 
representation 

In typical DFT+/7 approaches, the electronic struc¬ 
ture is given in terms of density matrices for each sub¬ 
space, e.g., the d shell. Unfortunately the orthogonal 
global axial representation which is most efficacious for 
computation is often not most convenient for analysis 
and understanding. This happens in calculations with 
non-trivial unit cells where inequivalent oxygen octahe- 
dra surround transition metal ions. Octahedral rotations 
and tilts mean the global axial system for the calcula¬ 
tion, here labelled x ', y ', z ', will differ from the native 
local axes, labelled x,y,z. Native axes for each octahe¬ 
dron point along the transition metal-0 bonds, and form 
the natural basis for understanding the electronic struc¬ 
ture of the transition metal d orbitals. We describe the 
details of a simple approach that rotates the density ma¬ 
trix, from the global to the local basis via polynomial 
transformations, with LaMnCU as our example. 

We choose a particular Mn ion and its nearest neigh¬ 
bor O atoms which identify an octahedral cage. Three 
Mn-0 bonds are chosen that point in approximately or¬ 
thogonal directions. The bonds are indexed i = 1,2,3, 
and we compute the difference vectors from the Mn to O 


defined by placing the unit vectors iq in the columns of R. 
It is at this point that we choose the ordering of the unit 
vectors to reflect the physical questions at hand. Note, a 
traditional choice is to align z with the apical bond, but 
other choices are possible: for example, in our work we 
have placed y along the non-Jahn-Teller ’apical’ Mn-0 
(see Figure [4]) while x and z span the Jahn-Teller active 
plane. 

This rotation represents a linear polynomial transfor¬ 
mation relating x',y',z' to x,y,z. The angular behav¬ 
ior of each d orbital is quadratic in the coordinates: 
3z' 2 — r' 2 , x' 2 —y' 2 , x'y', y'z', x'z', so it is straightforward 
to plug in and algebraically transform the polynomials to 
the unprimed (local) coordinate system. Performing the 
substitutions, using the orthogonal nature of the R ma¬ 
trix, and collecting terms, we find 


/ 3z 2 — r 2 \ 


/ 3z' 2 - r' 2 \ 

2 2 


2 n j2 

x — y 


x — y 

xy 

= c 

x'y 1 

yz 


y'z' 

V xz / 


\ x'z' / 


where the matrix C is 


C = 


( 2 (^ R 33 1) 

l( R2 31- R §2) 


6 R 31 R 32 
6R 3 2 R 33 
\ 6R 3 1 R 33 


H R i s~ r2 2 3 ) 

1 ( td2 _ jd2 _i_ td2 _ jd2 \ 

2 \- n 'll - rt 12 ' ■ rt 22 r *'2l) 

2(-Rn-Z?i 2 — R21R22) 
2(Rl 2 i?13 — R 22R23) 
2(RiiRi 3 — R21R23) 


\R13R23 

\{. R ll R 21 — R 12 R 22) 

R\1 R 22 + R12R2I 

R12R23 + R\3 R 22 
Rll R 23 + Rl3 R 21 


\ R 23 R 33 
{ R 21 R 31 ^ R 22 R 32) 
R 21 R 32 + R 22R 3 1 
R22R33 + R 23 R 32 
R 21 R 33 + R 23 R 31 


2 R 13 R 33 \ 

\ { R 11 R 31 — R 12 R 32) 

R11R32 + R12R3I 

R12R33 + Rl3 R 32 
R 11R33 + Rl3 R 31 / 


The matrix C is not unitary due to the fact that 
the bare polynomials 3z' 2 — r' 2 ,x' 2 — y' 2 ,x'y',y'z 1 ,x'z' 
are orthogonal but are not normalized. The normaliza¬ 
tion is done by averaging the squares of the functions 
(3z' 2 - r' 2 )/r' 2 ,(x' 2 - y' 2 )/r' 2 ,x'y'/r' 2 ,y'z'/r' 2 ,x'z'/r' 2 
over the surface of the unit sphere. We place these aver¬ 
ages, which are 4/5,4/15,1/15,1/15,1/15, respectively, 
on the diagonals of a diagonal scaling matrix S and 
then form the scaled and unitary transformation matrix 
D = S 1 / 2 C S _1/2 which is our final matrix relating the d 
orbitals in primed and unprimed coordinates. 


To give a feeling for how the method works, we take 
the experimental structure for LaMnOa crystal with a = 
5.736 A, b = 7.703 A and c = 5.540 A, as in Figure [I] 
In experimental structured LMO, consider the octahe¬ 
dron about the Mn atom at (0.00, 0.00, 2.77) A, which 
has two basal oxygens at O 1 = (1.12,-0.31,1.26) A 
and O 2 = (1.75,0.31,4.03) A, and an apical oxygen at 
O 3 = (—0.07,1.93,2.37) A. We form the normalized u; 
vectors, compute R and then C and upon normalization 
find 
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/ 0.06 

0.04 

0.01 

0.23 

0.31 


1 -0.03 

0.16 

0.88 

-0.24 

0.64 

D = 

0.58 

-0.74 

0.97 

-0.09 

0.27 


0.78 

0.54 

-0.71 

-0.19 

-0.33 


V -0.05 

-0.20 

-2.76 

-0.29 

0.99 


D can now be used to diagonalize the 5x5 density 
matrix in the sub-space of the Mn d orbitals. For a 
DFT+f7|J calculation with U = 8 eV and J = 2 eV, 
fixed at the experimental structure, the Mn d eigensys- 
tem is 



Here each eigenvector is a column vector with its eigen¬ 
value fie provided above it. Before rotation, it is hard 


to easily read off the nature of each eigenstate by inspec¬ 
tion. After rotation, the eigenvectors in the local basis 
are given by 



The local basis eigenvectors are clearly much “purer” 
as each vector has a component whose magnitude is 0.97 
or larger. And thus each configuration is easy to read 
off by inspection: the partially occupied state in the first 
column is essentially the x 2 — y 2 state while the last col¬ 
umn shows that the 3 z 2 — r 2 has become filled. We have 
strong orbital polarization in the e g manifold. This indi¬ 
cates the rotation to local octahedral coordinates success¬ 
fully diagonalized the eigensystem, and that the local ba¬ 
sis is physically relevant for understanding the electronic 
structure. 



